Adsorption-like Collapse of Diblock Copolymers 
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A linear copolymer made of two reciprocally attracting iV-monomer blocks collapses to a compact 
phase through a novel transition, whose exponents are determined with extensive MC simulations 
in two and three dimensions. In the former case, an identification with the statistical geometry of 
suitable percolation paths allows to predict that the number of contacts between the blocks grows 
like N 9 / 16 . In the compact phase the blocks are mixed and, in two dimensions, also zipped, in such 
a way to form a spiral, double chain structure. 
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Our knowledge of the conformational properties in so- 
lution of single, isolated polymers with inhomogeneous 
backbone sequence, is still extremely limited [Q. A ma- 
jor challenge in statistical physics consists in describing 
the possible conformational transitions of such systems, 
their universal scaling behaviors, and the nature of the 
different phases they connect, in finite dimensional situ- 
ations. In the last decades a similar program has been 
accomplished, to a substantial extent, for homopolymers 
while for heteropolymers, so far, most insight is lim- 
ited to approaches of mean field type g§. 

Block copolymers |^,|| are interesting as relatively ele- 
mentary members of the large family of heteropolymers, 
like proteins or polyampholytes, characterized by the 
above inhomogeneity. These copolymers are important 
as interface stabilizing agents, and can display intriguing 
phenomena of microphase separation, a major topic in 
soft condensed matter research. 

The status of the art, as far as studies of hetero- and co- 
polymer conformational transitions are concerned, sug- 
gests us to investigate in detail, without resort to mean 
field approximations, prototype problems involving rela- 
tively simple molecular architectures. Of particular inter- 
est will be those transitions which possibly reveal peculiar 
to systems with chain inhomogeneity, without counter- 
part in the homopolymer case. 

In this Letter we address one of such problems, the col- 
lapse from high temperature (T) swollen, to low T com- 
pact state of a copolymer whose two equally long blocks, 
A and B, attract each other with short-range forces. A 
physical realization could be that of oppositely charged 
A and B, immersed in a screening solvent. Another case, 
involving AA and BB attraction, rather than repulsion, 
is one in which the monomers of A are able to establish 
hydrogen bonds with those of B, which add to the Van 
der Waals forces, thus creating extra AB attractive in- 
teractions. In first instance we consider here the effect 
of AB attractive interactions acting alone. In both two 
dimensions (2D) and 3D, the collapsed phase has an ap- 



proximate alternating, periodic structure, as far as the 
space distribution of A and B monomers is concerned. 
We determine numerically, and predict exactly in 2D, 
the exponents of this collapse, showing that it has qual- 
itative features in common with adsorption phenomena 
Q. Nevertheless, the universality class coincides with 
neither that of adsorption [0 , nor that of theta collapse 
of homopolymers 6 . 

We model the diblock copolymer by a self-avoiding 
walk (SAW) w of length \w\ = 2N steps (monomers) 
on square and cubic lattice. The walk w consists of N 
consecutive monomers of type A (wa), followed by N 
monomers of type B (wb)- The Hamiltonian takes the 
form 



H{w) = - 



(1) 



where e is a positive energy, and the sum is supposed 
to run over pairs < i, j > of lattice sites visited by the 
copolymer ( the AB junction excluded), which are near- 
est neighbors (n.n.). We will also consider the possible 
inclusion in H of extra terms representing analogous n.n. 
AA and BB repulsive or attractive potentials R. 

The scaling regime of a polymer can be conveniently 
described by the v exponent governing the asymptotic 
behavior of the canonical average radius of gyration, 



(R g ) 



T,w exp(-H(w)/T)R(w) 
J2 w exp(~H(w)/T) 



(2) 



where the sums extend to all 27V-step configurations w 
of the copolymer, with radius R(w) relative to the cen- 
ter of mass. In the swollen, high T regime one should 
expect v = vsAWi with vsaw = 3/4 and 0.588 (ljj 
in 2D and 3D, respectively. In a compact, low T phase, 
v = l/d. Another important physical quantity is the 
specific heat C — -kd < H > /dT, which at a confor- 
mational transition between high T and low T regimes is 
expected to obey a scaling of the form: 



C ~ N^-^FHT -T C )N^) 



(3) 
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for large N and for T close to the transition temperature 
T c , F being a suitable scaling function. For T — T c , v 
should take a peculiar value, v c . 

By extensive Monte Carlo sampling we computed (R g ), 
C and (H) of the copolymer for N up to 800 in both 2D 
and 3D. The sampling was based on a multiple Markov 
chain (MMC) method, which has been shown to be par- 
ticularly suitable to deal with polymers at low T ||. 
Several Markov chains, each one designed to sample at 
a different T, are generated in parallel, by using both 
pivot H and local moves p0[ . The sampling at low T 
is then considerably enriched by swapping configurations 
between chains contiguous in T. Since each Markov chain 
is ergodic, so is the composite one J§| . The swapping pro- 
cedure dramatically reduces the correlations within each 
chain, and produces little CPU waste time since, in any 
case, one is interested in obtaining data at several T's [[§. 
Both in 2D and 3D the simulations have been performed 
by using 25 parallel chains covering a range from T = oo 
to deeply inside the collapsed phase (T = 0.5). 

Our data for C ( Fig.Q) signal clearly the presence of 
a transition, with peaks in the T-dependence sharpen- 
ing and growing with N, consistently with Eq.(3). The 
corresponding T c and the crossover exponent <f> could be 
deduced from the N dependence of the height, h(N), 
and position, T C (N), of these peaks. Indeed, in the scal- 
ing limit we expect h(N) ~ TV 2 *" 1 and T C (N) -T c (oo) ~ 
N~^, for increasing N (T c (oo) = T c ). 




FIG. 1. Specific heat as a function of T for TV = 50 (*), 100 
(❖), 200 (d),300 (a), 400 (•), 600 (o) and 800 (filled triangles). 

Since a linear least squares fit of the log of the asymp- 
totic form of h gives a very large % 2 statistical error, 
we consider h(N) = AiV 2 * _1 (l + B/N) where a scal- 
ing correction 1/N is included. The least squares fit in 
this case gives a lower x 2 and we obtained (see FigJ^) 
(j) = 0.56 ± 0.02 in 2D and cf> = 0.60 ± 0.01 in 3D. These 
values of cf> allowed to extrapolate T c (oo)/e = 1.5 ± 0.2 
and T c (oo)/e = 2.9 ± 0.2, in 2D and 3D, respectively. 
From log-log plots of (R g ) vs. N we also estimated the 
v exponents at several temperatures. At T — T c , as de- 
termined above, we obtain v c = 0.748 ± 0.008 in 2D, 
and v c — 0.583 ± 0.007 in 3D. In fact v stays approx- 



imately constant for T > T c . Only for T < T c , v de- 
creases, approaching values consistent with a compact 
phase [y ~ 1/d) for low T. 
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FIG. 2. Fits of h(N) in 2D (triangles) and 3D (squares). 

The v c estimates appear fully compatible with the ex- 
ponents appropiate to a swollen SAW. In this respect 
the new collapse is analogous to a polymer adsorption 
on an attractive, impenetrable wall ||. Indeed, at an 
adsorption transition one also finds v a ds = ^saw- How- 
ever, the crossover exponents 4> determined here are def- 
initely different from those describing the growth with 
N of the number of polymer-wall contacts in adsorption 
{(j) ads = 1/2 and (j) ads = 0.496 ± 0.004 in 2D § and 3D 
pT[ , respectively) . Further analogies with adsorption are 
revealed by the large N behavior of the average energy 
(or (Nab)) per monomer. FigJ| indicates crossings of the 
various (H)/N curves. The temperature range at which 
the various crossings concentrate is a signal of the tran- 
sition and is consistent with the above estimate of T c . 
Moreover (H) /N -> for T > T c and (H) /N -> const 
for T < T c (inset of Fig. ||). This behavior is completely 
different from what is found for the theta collapse, where 
the energy density curves do not cross each other, and 
(H)/N is asymptotically nonzero for all temperatures. 
The average number of polymer-wall contacts at oppo- 
site sides of an adsorption transition behaves similarly to 
(Nab) here. In addition, like in adsorption, at T = T c 
we find that (N AB )/N = -(H) /N - N^ 1 . By fit- 
ting this behavior, we obtained alternative 4> estimates 
(4> = 0.55 ± 0.04 and = 0.61 ± 0.03, in 2D and 3D 
respectively), less sharp, but fully consistent with those 
based on the specific heat. 

Apart from being a transition from swollen to compact 
state, the diblock copolymer collapse has nothing in com- 
mon with a homopolymer theta collapse, which is char- 
acterized by a peculiar v exponent different from that of 
a high T swollen chain (v® — 4/7 and 1/2, in 2D and 3D, 
respectively |||^]). Also the theta point crossover expo- 
nent is clearly different from the cf>'s here (cf>Q = 3/7 in 2D 
H and 4>q = 1/2 in 3D M). The copolymer collapse is 
indeed a novel transition, determined by the peculiar in- 
homogeneous backbone structure and by its interactions. 
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This and other remarkable features discussed below, con- 
fer this transition a prototypical interest, also beyond its 
possible applications. 




FIG. 3. (H)/N vs T in 3D for N = 50 (*), 100 (o), 200 
(□),300 (a), 400 (•), 600 (o) and 800 (filled triangles). The 
inset shows the convergence of (H) /N to for T > T c and to 
a constant for T < T c . 




FIG. 4. Typical configuration at low T (T — 0.5 < T c ) in 
d=2 for N = 100. Filled and empty circles denote, respec- 
tively, A and B monomers. 

Further insight can be gained by considering the struc- 
ture of the dense phase. In both 2D and 3D, typical dense 
configurations generated at low T by our sampling tech- 
nique present a sort of maximal alternance of A and B 
monomers, compatibly with the constraint of forming a 
diblock chain. In 2D (Fig. [|) A and B are also zipped 
on each other, and the resulting double chain structure 
forms spirals, making the number of attractive interac- 
tions maximal. The zipped, double chain represents a 
concrete realization of a polymer with orientation depen- 
dent interactions. For models of such polymers, spiral 
compact phases have been recently observed in 2D |Q. 

In 3D the compact phase has more complex configura- 
tions and substantial entropy. In this case zipping into a 
double chain does not occur. On the other hand, the A 
and B monomers are not segregated, and their alternance 
in space resembles, to some extent, an antifcrromagnetic 
microphase structure. By assuming v c — 3/4, we can 
predict 4> exactly in 2D, following a strategy of identifi- 
cation with geometric percolative objects, which in the 
past revealed extremely fruitful for homopolymer prob- 
lems P jT^JT^ ]. The idea is that the statistics of suitably 



chosen percolative contours could reproduce the critical 
behavior of interacting polymers in the conditions of in- 
terest. The homopolymer at the theta point in 2D is 
known to have the same statistics as the external perime- 
ter (hull) of a critical percolation cluster p7| , and this 
coincidence is at the basis of the full exact characteri- 
zation of theta scaling j|,[l^,[ll| . O n square lattice, the 
incipient infinite cluster of the percolation problem for el- 
ementary squares ( only squares with an edge in common 
are connected) is expected to have an externally acces- 
sible hull which, besides assuming the configurations of 
a self-avoiding ring, has the statistical fractal dimension 
Dd = 1/vsaw = 4/3 of a swollen SAW @. 

This has been recently established on a rigorous basis 
in Ref. JllJ . By externally accessible hull we intend that 
there should exist at least one path of connected vacant 
squares by which any neighborhood of the perimeter can 
be connected to points at infinite distance on the lattice. 
We try to identify our A and B blocks at the transition 
with the two halves in which the externally accessible hull 
of the cluster is divided by two diametrally opposite, dis- 
tant points (Fig. ||). The identification of A and B with 
the copolymer blocks at collapse will make sense provided 
the reciprocal n.n. contacts of these halves have a posi- 
tive fractal dimension, matching the correct <p exponent, 
as discussed below. This positive fractal dimension is in 
fact the result of an effective attractive AB interaction. 
As sketched in Fig. ||, the fact that A and B are identi- 
fied with externally accessible hull portions excludes from 
the count the full hull self-contacts possibly produced by 
the presence of reentrances inaccessible from the exterior. 
Indeed, such reentrances can not be included in A or B, 
and it is in virtue of their exclusion that the accessible 
perimeter has a fractal dimension D c \ < Dhuii = 7/4. 
In Ref. jl9| Coulomb gas results for the T = loop gas 
are used for a systematic determination of the dimen- 
sions of various fractal sets in the cluster of the hexagon 
percolation problem. For this problem it turns out to be 
more easy to discuss geometry and effective interactions 
of these sets. In view of the peculiar connectivity proper- 
ties of the hexagonal lattice, hull reentrances are always 
accessible by at least one path of vacant hexagons. Thus, 
external accessibility requires now existence of at least 
two distinct vacant paths joning the neighborhood of the 
hull to infinity [HJ. The probability with respect to all 
percolative configurations of a given hull profile, can be 
interpreted as the result of a decision process (between 
being occupied or vacant) for the hexagons which are 
progressively touched by the hull itself. It is thus easy to 
realize that AB n.n. contacts have the same effect as at- 
tractive interactions in a polymer problem. Indeed, while 
normally each step of the hull involves decision for a new 
hexagon (and thus a probability factor 1/2, the hexagon 
percolation threshold, in the total weight), at an AB con- 
tact the involved hexagon has already made his decision, 
and no factor is required. Thus, AB contacts increase the 
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probability of the hull configuration. Besides determin- 
ing D c i exactly, the methods of Ref. [jl9| allow us to con- 
clude that the set of AB contacts described above must 
have a fractal dimension D c2 — 3/4, which coincides with 
that of the "red" hexagons (or squares in our case) of the 
cluster, i.e. of those hexagons whose removal would inter- 
rupt the cluster connection between the opposite distant 
points |^0| . D c i and D c2 are also expected to be universal 
with respect to different lattice structures (ll]] , and thus 
should hold also for square lattice. The proposed identi- 
fication leads to 4> = 9/16: indeed, for iV-steps accessible 
half-perimeters, the average number of contacts should 
scale like N D - 2 / D ^. D c2 /D cl = 9/16 = 0.56.. is strik- 
ingly close to our numerical determination, supporting 
the correctness of our assumptions [^lj . 




FIG. 5. Points 1 and 2 separate the accessible hull portions 
A and B. The light- continous line represents an inaccessible 
reentrance of the hull, k : red square corresponding to AB 
n.n. contacts, m : contacts of A with a hull reentrance. The 
union of A and B has fractal dimension 4/3. By including in 
the hull definiton also inaccessible reentrances, the dimension 
becomes 7/4. 



In summary, we gave a description of the collapse tran- 
sition undergone by a copolymer with two reciprocally 
attracting blocks. We could determine with high preci- 
sion, and, by an identification with percolation paths, ex- 
actly in 2D, the universal exponents of the new transition, 
which has analogies with homopolymer adsorption and 
theta collapse, but falls in a universality class different 
from both. The low T dense phase has antifcrromagnetic 
features and, in 2D, turns out to be spiral, with zipped 
A and B. To our knowledge, the exact determination of 
<j) in 2D is the first application of percolation results to a 
genuinely inhomogeneous polymer problem, and demon- 
states the potential of methods which, about a decade 
ago, led to the correct characterization of the homopoly- 
mer theta point. We checked that the universality class 
of the transition remains the same if, more consistently 
with the physical picture of oppositely charged blocks, 
AA and BB repulsions are included in Eq.(l). A natural 
extension of our study consists of a systematic explo- 
ration of the phase diagram in regions where, e.g., the 
AA and BB interactions are also attractive and compete 
with the AB one Q| . This corresponds to a copolymer in 
which the difference between AA ( or BB) and AB inter- 



actions is produced by hydrogen bonds. We preliminarily 
verified that this competition generates a very rich dia- 
gram, with both segregated and unsegregated compact 
phases and interesting new transitions. 
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